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(54) Title: METHOD FOR PREDICTING DE NOVO BIOMACROMOLECULE CRYSTALLIZATION CONDITIONS AND FOR 
CRYSTALLIZATION OF THE SAME 

(57) Abstract: A method de novo is provided for predicting crystallization conditions and for crystallizing biomacromolecules, in 
particular proteins. The method provides a simple, quick and precise approach in determining the biomacro molecule solubility in 
different solutions, as well as the boundary between crystallization and aggregation. Because the method relies only on monitoring 
the assembly behavior of the biomacro molecule at the surface of a solution, it has general applicability and requires a relatively 
short amount of time to provide results that are reliable. Because there is no need to first crystallize the biomacromolecule, smaller 
amounts of protein suffice as compared with amounts required for crystallization. Because the method works by measuring the 
surface tension or surface pressure of the surface of the biomacromolecule solution, it is easy, precise and quick. Furthermore it is 
cost-effective in requiring simple and inexpensive equipment. Additionally, the method not only provides information on the solvent 
conditions where crystallization, rather than amorphous aggregation, takes place, but also on the critical equilibrium condition of 
the protein. This enables the user to restrict experimental parameters to protein concentrations above the equilibrium value, that are 
viable for crystallization. 
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METHOD FOR PREDICTING DE NOVO BIOMACROMOLECULE CRYSTALLIZATION 
CONDITIONS AND FOR CRYSTALLIZATION OF THE SAME 

Field of Invention 

This invention relates to a de novo method of predicting crystallization conditions. In 
particular, the invention relates to a method of predicting the conditions for crystallizing 
biomacromolecules from solution, and more specifically for crystallizing proteins from 
solution. 

Background 

In most cases of protein crystallization experiments, the final product is not a single 
crystal but amorphous aggregation. To predict the likelihood of the formation of either a 
crystal or an amorphous aggregation, the second virial coefficient B22 is customarily employed 
by many groups [1 — 5]. Serving as an indicator of intermolecular interactions, ,822 is positive 
when these interactions are repulsive, and negative when these interactions are attractive. As a 
consequence, a necessary condition for crystallization is obtained when the second virial 
coefficient B22 lies in a so-called "crystallization window": - 8 xlO" 4 < B22 < - 2 x 10" 4 ml 
mol/g 2 . 

The criterion based on the second virial coefficient has its advantages, as it gives a 
discriminating response. However this criterion does not work in some cases in which B22 
could well lie within the crystallization window, but the experiment gives amorphous 
aggregation. One of the reasons for this failure is that the criterion provided by the second 
virial coefficient takes into account only the interactions between the biomacromolecules. 
However, such intermolecular interactions determine biomacromolecule crystallization only 
partially. Apart from this, the B22 criterion only attempts to determine the solvent conditions 
for which crystallization rather than amorphous aggregation would take place. However, for 
biomacromolecule crystallization to be possible in the first place, an additional condition must 
be satisfied, namely that the biomacromolecule concentration must exceed its equilibrium 
value in the context of the applied solvent conditions. The B22 criterion does not address this 
problem, and gives no information whatsoever on the critical equilibrium condition of the 
biomacromolecule. 
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The crystallization of biomacromolecules involves a nucleation and growth process, 
determined to a large extent by kinetics. Kinetics refers to the way biomacromolecules move 
in a solution, the rate at which they are transported, and the way they are incorporated in the 
biomacromolecule crystals at the crystal surface. The crystallization window provided by the 
second virial coefficient disregards kinetic and other factors, that are unrelated to 
intennolecular interactions but nevertheless largely influence crystallization. 

Alternative methods have been developed to circumvent the drawbacks presented by 
the use of the second virial coefficient in attempting to predict crystallization. These methods, 
the so-called high throughput screening methods, are solely empirical involving a large 
number of solution matrices. These methods are costly, cumbersome, and time-consuming; 
they require large investments in expensive robots. As protein crystallization is normally a 
lengthy process and can be affected by about 20 different parameters, these methods have in 
many cases no general applicability and suffer from a low success rate. 

What is needed is a de novo method offering a twofold advantage above the prior art. 
Firstly, such a method should determine the equilibrium biomacromolecule concentration, so 
that the values of the biomacromolecule concentration used in the process of crystallization in 
the prevailing experimental conditions can be restricted to values above the equilibrium value. 
Secondly, such a method should predict more reliably the crystallization conditions referring 
to the solvent experimental parameters. Finally, such a method should have a more general 
applicability and be simple, easy to apply, not wasteful on the biomacromolecule or the 
protein. It should not require complex equipment, or rely on heavy investments for its 
application, and it should be amenable to being utilized by both the institutional and industrial 
establishments, while at the same time offering substantial automation advantages, and being 
much less costly than the methods recited in the prior art. Such a de novo method does not 
exist in the prior art. 
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Summary 

Ade novo method for predicting crystallization conditions for biomacromolecules is 
provided, giving a more reliable prediction criterion than has been possible thus far. Unlike the 
second virial coefficient B22 used to predict crystallization conditions based solely on 
intermolecular interactions, the present invention combines information on both the 
intermolecular interactions and kinetic effects to prescribe crystallization conditions. The 
improvement above the empirical methods available in the prior art and the prediction 
methods provided by calculating the second virial coefficient, is due to the incorporation of the 
kinetic effects that largely determine biomacromolecule crystallization, and offers the 
following advantages: (1) it establishes an equilibrium macromolecule concentration that 
should be exceeded by the applied macromolecule concentration used for crystallization; (2) 
only small amounts of protein are required; (3) the measurements are much quicker than either 
through use of the B22 approach or through use of the empirical screening methods; (4) the 
simplicity of the required measurement of the assembly parameter (e.g. surface tension) of the 
biomacromolecule solution facilitates the handling in an ordinary laboratory. 

One object of the present invention is to identify the phase boundary between a crystal 
phase and a liquid phase by determining the biomacromolecule solubility, under prevailing 
experimental conditions. Another object of the invention is to identify the boundary between a 
crystallization and an aggregation under prevailing experimental conditions. 

Therefore in accordance with a first aspect of the invention there is disclosed: 
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A method for predicting a crystal equilibrium condition for biomacromolecule crystallization 
and for crystallizing a biomacromolecule, comprising 

setting up at least one biomacromolecule solubility experiment comprising the steps of 

a) preparing a solution of the biomacromolecule in a solvent, the solution having a 
biomacromolecule concentration, 

b) selecting a variable quantity, 

c) selecting an assembly parameter, 

d) monitoring a response of the assembly parameter while varying the variable quantity in 
a suitable way so that the response exhibits a transition, 

e) obtaining an equilibrium biomacromolecule concentration based on the transition, 

f) defining a crystal equilibrium condition according to which a biomacromolecule 
crystallization concentration exceeds the equilibrium biomacromolecule concentration, 

and crystallizing the biomacromolecule. 

In accordance with a second aspect of the invention there is disclosed: 
A method for predicting an aggregation boundary condition for biomacromolecule 
crystallization and for crystallizing a biomacromolecule, comprising 
setting up at least one aggregation boundary condition experiment comprising 

a) preparing a solution of the biomacromolecule, 

b) selecting a variable quantity, 

c) selecting an assembly parameter, 

d) measuring the assembly parameter at different times, 

e) registering an equilibrium assembly parameter 

f) deriving a crystallization coefficient from the equilibrium assembly parameter, the 
crystallization coefficient being associated with the variable quantity, 

g) using an aggregation indicator to define an aggregation boundary condition for the 
biomacromolecule, the aggregation boundary condition prescribing that an aggregation 
occurs when the crystallization coefficient associated with the variable quantity is larger than 
the aggregation indicator, 

and crystallizing the biomacromolecule. 

Brief Description of the Drawings 

Embodiments of the present invention are illustrated by the following drawings: 
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FIG. 1 is an illustration of predictive conditions; 

FIG. 1A is an illustration of a crystal equilibrium condition, leading to an equilibrium 
biomacromoiecule concentration; 

FIG. IB is an illustration of an aggregation boundary condition, leading to a crystallization 
coefficient; 

FIG. 2 is an illustration of a crystallization process from a solution of biomacromolecules, 
involving a transition from a liquid phase to an orderly solid phase, that is a crystal; 
FIG. 3 is an illustration of an amorphous aggregation process involving the transition from the 
liquid phase to a disorderly solid phase; 

FIG. 4 is an illustration of analogous processes to the crystallization and the amorphous 
aggregation that take place in the solution; 

FIG. 5 is an illustration of a response of an assembly parameter to a tendency of the 
biomacromolecules to assembly; 

FIG. 6 is an illustration of the way a solubility curve relates to the crystal equilibrium 
condition. 

FIG. 7 shows plots of the assembly parameter against time, when the solution of 

biomacromolecules includes a protein lysozyme and a salt sodium chloride; 

FIG. 8 shows plots of the assembly parameter against the square root of the time; 

FIG. 9 shows plots of the logarithm of the time derivative of the assembly parameter; 

FIG. 10 is an illustration of a time-dependence profile of a macromolecule behavior; 

FIG. 1 1 shows plots of a diffusion time, a penetration time and a rearrangement time; 

FIG. 12 shows a histogram illustrating how crystallization coefficients relate to the 

aggregation boundary condition by prescribing an aggregation indicator; 

FIG. 13 shows a plot of the assembly parameter against time, when the solution of 

biomacromolecules includes the protein lysozyme and a salt ammonium sulfate; 

FIG. 14 shows a plot of the assembly parameter against time, when the solution of 

biomacromolecules includes a protein concalavine A and the salt ammonium sulfate; and 

FIG. 15 shows plots of the assembly parameter against time, when the solution of 

biomacromolecules includes a protein bovine serum albumin (BSA) and the salt ammonium 

sulfate. 

Detailed Description 
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The primary features of a first embodiment of the invention of a method of predicting 
biomacromolecule crystallization conditions and for crystallizing biomacromolecules are 
provided hereinafter with reference to FIG. 1A in FIGs. 2-6. The primary features of a second 
embodiment of the invention are provided hereinafter with reference to FIG. IB in FIGs. 2-5 
and FIGs. 7-15. In the first embodiment, a crystal equilibrium condition in FIG. 1A is 
expressed by means of a macromolecule solubility curve serving as a boundary separating two 
regions of experimental parameter values: a region where crystallization can occur, and a 
region where crystallization cannot occur. 

The method described in the first embodiment establishes a biomacromolecule 
equilibrium concentration in the context of the applied experimental conditions. The 
biomacromolecule concentration to be used for crystallization must exceed the obtained 
equilibrium value. In a second embodiment, an aggregation boundary condition in FIG. IB is 
expressed by means of a window of experimental parameter values above which the 
amorphous aggregation is likely to occur. The terms "aggregation" and "amorphous 
aggregation" are used interchangeably. 

Theoretical Introduction 

In attempts to crystallize biomacromolecules from a solution, it is desirable to obtain as 
much as possible single crystals with as few defects as possible, and to avoid amorphous 
aggregations of molecules, since amorphous aggregates are not crystals. The present invention 
takes advantage of the property of biomacromolecules to have mixed hydrophobic and 
hydrophilic regions. This property results in a tendency for these molecules to assembly either 
in the bulk or at the surface of the solution. In this disclosure, the surface of the solution can be 
adjacent to another material or to empty space, and hence the surface can be in contact with a 
solid or with a liquid or with a gas, that is usually air. The surface of the solution has a surface 
tension and a surface pressure, which terms in this case include an interfacial tension or an 
interfacial pressure. 

It is possible to define one or more assembly parameters that reach a critical response 
as increasingly more molecules participate in assembly formation. For example, the tendency 
of biomacromolecules to assembly in a solution can be monitored by taking density, 
conductivity, detergency, osmotic pressure, surface tension or surface pressure measurements 
of the solution. 



WO 2005/080422 PCT/SG2005/000051 

-7- 

Biomacromolecule crystallization conditions are reflected in the tendency of 
biomacromolecules to assembly. The crystallization conditions and the assembly 
characteristics of biomacromolecules are governed by both the intermolecular forces and the 
kinetic effects [6-8], depending on the experimental situation. The present invention prescribes 
two procedures of measuring the assembly parameters in solution so as to determine reliable 
crystallization conditions without the prior need to carry out crystallization experiments. 
Without limitation to the scope of the invention, the present examples illustrating the 
invention show that the disclosed method is particularly applicable when the 
biomacromolecule is not prone to severe unfolding at the surface of the solution or at an 
air/solution interface. 

Application 

In FIG. 2 a crystallization process from the solution 120 comprising 
biomacromolecules 130 involves a transition 122 from a liquid phase 124 to an orderly solid 
phase, that is a crystal 126, wherein solid-phase biomacromolecules 132 are regularly arrayed 
and regularly oriented. The biomacromolecules 130 in the liquid phase 124, become adsorbed 
biomacromolecules 136 on a surface 138 of the crystal 126, and subsequently they acquire a 
suitable orientation 140 in order to become incorporated as in the crystal 126 as the solid- 
phase biomacromolecules 132. The crystallization occurs because a transport time, being an 
amount of time required for the biomacromolecules 130 in the solution 120 to reach the 
surface 138 of the crystal 126, otherwise known as a diffusion time r d ifr, is longer than an 
integration time t mt9 wherein the integration time is an amount of time required for the 
adsorbed biomacromolecules 136 to become the solid-phase biomacromolecules 132, and it 
includes an amount of time required for their rearrangement. 

In FIG. 3 on the other hand, the amorphous aggregation involves the transition 122 
from the liquid phase 124 to a disorderly solid phase 128, wherein aggregate 
biomacromolecules 134 have irregular positions and irregular orientations. The amorphous 
aggregation occurs because the diffusion time r d ifr required for the biomacromolecules 130 in 
the solution 120 to reach the solid phase 128 is shorter than the time required for their 
rearrangement and integration. 

In FIGs. 2 and 3 the solution 120 has an appropriate temperature and an appropriate 
pH. In the invention there is no restriction on the appropriate temperature, which is usually 
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taken to be around a room temperature. A buffering material may be introduced in the solution 
120 to maintain a desirable level of the pH. The solution 120 can also comprise additives like 
small or large organic molecules or salt. By way of example, the biomacromolecule 130 in 
FIGs. 2 and 3 is a protein, more specifically lysozyme (about 14.3 kDalton) or concalavine A 
(about 23 kDalton), and the additive is a salt, more specifically sodium chloride or ammonium 
sulfate. 

FIG. 4 illustrates that processes analogous to the crystallization or the amorphous 
aggregation take place in the solution 120 containing the biomacromolecules 130. The solution 
120 is placed in a container 140, and the solution has a surface 152, which in this application 
is adjacent to air. During a diffusion step 142, the biomacromolecules 130 take an average 
amount of time, that is the diffusion time / diff , to diffuse 148 towards the surface 152. 
Subsequently, during a penetration step 144, the biomacromolecules 130 take an average 
amount of time, that is a penetration time to penetrate 150 the surface 152. Subsequently, 
during a rearrangement step 146, the biomacromolecules 130, that penetrated 150 the surface 
152 take an average amount of time, that is a rearrangement time t m , to undergo a 
configuration^ rearrangement 154 and finally become integrated in the surface 152. Thus the 
integration time t mt is the sum of the penetration time and the rearrangement time, t mt = t pm + 
/arr, because it equals an amount of time required to complete the penetration step 144 and the 
rearrangement step 146. 

In FIG. 5 the biomacromolecules 130 have a biomacromolecule concentration 156, and 
the additives have an additive concentration 157. Without limitation to the scope of the 
invention, in the experimental runs of the first embodiment of the invention, the 
biomacromolecule concentration 156 is taken, e.g. 4 mg/ml lysozyme, while the additive (e.g. 
salt) concentration 157 is allowed to vary. In some other experimental runs of the first 
embodiment of the invention, the additive concentration 157 is taken, e.g. 1 M (as in FIG. 5) 
or 2 M salt, while the biomacromolecule concentration 1 56 is allowed to vary. 

Referring to FIGs. 5 an 6, in the first embodiment of the invention, measurements of an 
assembly parameter 1 64 (for example a surface tension or a surface pressure) with respect to a 
variable quantity 158 (for example the biomacromolecule concentration 156, the additive 
concentration 1 57, the pH or the temperature), are utilized to express the crystal equilibrium 
condition. According to that condition, as will be explained further on, in FIG. 1 A (FIG. 6), 
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the crystallization is expected in a region 172 delineated by (and lying above and to the right 
of) a biomacromolecule solubility curve 170. 

In the second embodiment of the invention, a probability that the aggregation (FIG. 3) 
will dominate above the crystallization (FIG. 2), increases when the integration time t mt is 
substantially longer than the diffusion time f d iff. That probability can be quantified as 
illustrated in FIG. IB (FIG. 12), by means of a crystallization coefficient 232 ft^t = /int/*diff 
defined as a ratio of the integration time to the diffusion time. Good quality crystallization 
occurs when fe^t is in a suitable range. The higher the value of Abiyst, the higher the probability 
that the amorphous aggregation will occur. When above A«ys t <1, the solution 120 containing 
the biomacromolecule 130 exists in a so-called "steady state", in which neither the 
crystallization (126 in FIG. 2) nor the aggregation (128 in FIG. 3) can take place. The 
aggregation indicator resulting from the crystallization coefficient should therefore certainly 
lie well above —1; the likely reason for this is that if A^iyst <1, the diffusion step 142 of the 
biomacromolecules 130 in the solution 120 lasts such a long time, that it becomes difficult for 
the biomacromoles 130 to come in contact and group with one another. In practice however, in 
our experimental applications, it turns out that a realistic minimum value for fccryst, allowing 
biomacromoles 130 to come in contact and group lies around 4 or 4.5. 

Because the crystallization coefficient 232 is a dimensionless ratio, it is expected to 
have general applicability for most biomacromolecules. The method is illustrated in this 
experiment by choosing as examples the protein lysozyme to serve as a model system for the 
biomacromolecule crystallization, as well as the protein concanavalin A. It is therefore 
expected that the aggregation indicator obtained from the crystallization coefficient derived for 
the proteins lysozyme and concanavalin A will serve as a standard criterion to define the 
aggregation boundary condition for the crystallization of most biomacromolecules. 

In the second embodiment of the invention, the solution 120 in FIGs. 4 and 5 is fairly 
dilute. Without limitation to the scope of the invention, in some experimental runs the 
biomacromolecule concentration is taken 156, 158 in FIGs 4 and 5 is taken in the range 0.01 
mg/ml to about 1.2 mg/ml. Measurements of the diffusion time / diff , the penetration time t pfm 
and the rearrangement time are used to calculate the crystallization coefficient 232 kc^t = 
Wfaff which in turn is used to express the aggregation boundary condition, as explained 
further below with reference to FIGs. 7—15. 
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First Embodiment: Crystal Equilibrium Condition 

In the first embodiment of the invention, in FIGs. 4 — 6, the crystal equilibrium 
condition is determined by constructing the biomacromolecule solubility curve 170. 
FIG. 5 illustrates a response of the assembly parameter 164 to the activity of the 
biomacromolecules 130 in the solution 120. A suitable assembly parameter 164 for this 
application is the surface tension or the surface pressure, and in this case it is taken to be the 
surface tension, that can be measured with a tensiometer. The surface pressure p is related to 
the surface tension s by the equation: p — s-sq 9 where s 0 is the surface tension of the solvent 
in the absence of biomacromolecules. A suitable variable quantity 158 for this application is 
the biomacromolecule concentration 156 or the additive concentration 157, the pH, or the 
temperature. When a suitable variable quantity 158 has been selected, it can be made to vary in 
order to sample the response of the assembly parameter 164. In this case the variable quantity 
158 is taken as the biomacromolecule concentration 156, whereas the additive concentration 
157 is taken to be 1 M salt, e.g. NaCL 

FIG. 5 shows a plot of the assembly parameter 164 (the surface tension) against the 
logarithm 159 of the variable quantity 158 (the biomacromolecule concentration 156), 
according to which the assembly parameter 164 decreases as the variable quantity 158 is 
gradually increased. This is illustrated by the difiusion step 142, followed by the penetration 
step 144, subsequently resulting in completion of the rearrangement step 146. Upon 
completion of the rearrangement step 146, a complete layer of rearranged biomacromolecules 
has been assembled at the surface 152, allowing no space for any more biomacromolecules 
130 to penetrate the surface 152. Therefore a further increase of the variable quantity 158, 
being in this experiment the biomacromolecule concentration, cannot cause a further decrease 
in the assembly parameter 164, being in this experiment the surface tension, but will instead 
result in accumulation of biomacromolecules 130 in the bulk of the solution 120. 

FIG. 5 illustrates a biomacromolecule solubility experiment. Thus in FIG. 5 as more 
and more biomacromolecules 130 tend towards assembly formation, the response of the 
assembly parameter 164 shows a transition 162 between a changing response 161 and a 
substantially constant response 163 of the assembly parameter 164. The transition 162 is 
associated not only with a critical response of the assembly parameter 164, but also with a 
critical magnitude 165 of the variable quantity 158. In this case the critical response of the 
assembly parameter 164 is equal to the surface tension of 53 mN/m. The critical magnitude 
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165 of the variable quantity 158 is equal to the lysozyme concentration of 4 mg/ml; this 
corresponds to the logarithm 159 of the variable quantity 158 equal to 1.4, when the additive 
concentration 157 is equal to 1 M NaCl. 

Because in this particular example the assembly parameter 164 is the surface tension, it 
decreases as the variable quantity 158, that is the biomacromolecule concentration 156, is 
increased, and hence the critical response of the assembly parameter 164 is substantially 
minimal. The crystal equilibrium condition can be expressed by means of the critical 
magnitude 165 of the variable quantity 158: when the temperature and pH are held at their 
predetermined values, and when the salt concentration is held at 1 M NaCl, then the 
crystallization cannot occur for values of the lysozyme concentration 156 falling below the 
equilibrium biomacromolecule concentration, the value of which in the prevailing 
experimental conditions is 4 mg/ml. A similar critical behavior is observed when the additive 
concentration 157 is allowed to vary, while the biomacromolecule concentration 156 is held 
constant at 4 mg/ml lysozyme (drawing not shown). 

The aforementioned behavior should not be construed to be typical of a general case 
covered by the scope of the invention in which the transition 162 refers to the changing 
response and the substantially stable or the substantially unchanging response, considering that 
the changing response need not imply a purely decreasing or a purely increasing response. 

FIG. 6 illustrates how the biomacromolecule solubility curve 170 can be constructed in 
order to express the crystal equilibrium condition. The solubility curve 1 70 follows when two 
of the following quantities: the biomacromolecule concentration 156, the additive 
concentration 157, the pH and the temperature, are allowed to form a pair of the variable 
quantities. In the example of FIG. 6 the biomacromolecule concentration 156 and the additive 
concentration 157 form the pair of the variable quantities, while the pH is kept fixed. A series 
of transitions 1 62 is obtained, each transition associated with the critical response of the 
assembly parameter 164. In the particular experiment described here, six transitions 162 were 
measured and employed in order to construct the solubility curve 170. For example, a point 
171 on the solubility curve 170 indicates the following pair of the variable quantities: the 
biomacromolecule concentration 156 is 4 mg/ml lysozyme while the additive concentration 
157 is around 1 M NaCl. 

The corresponding crystal equilibrium condition follows by specifying that the 
crystallization can occur in a crystallization region 172 in which either one of the pair of the 
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variable quantities assumes values at or above the corresponding critical magnitude 165 on the 
solubility curve 170. Therefore the crystal equilibrium condition as illustrated in FIG. 6 is 
obtained by recording a series of the critical magnitudes 165 of the variable quantity and the 
corresponding critical response of the assembly parameter 164, as illustrated in FIG. 5. 

When the biomacromolecule concentration 156 is used as a variable quantity, the 
critical magnitude 165 is the equilibrium biomacromolecule concentration. When some other 
parameter is used as a variable quantity, e.g. one of the additive (salt) concentration 156 or the 
pH, or the temperature, the equilibrium biomacromolecule concentration is, in the context of 
the employed experimental conditions, given by the applied biomacromolecule concentration 
1 56. In either case the crystal equilibrium condition prescribes that for crystallization to take 
place, the crystallization biomacromolecule concentration must exceed the equilibrium 
biomacromolecule concentration resulting from the biomacromolecule solubility experiment. 

For the measurement of the points on the solubility curve the Wilhelmy plate method is 
employed using a K14 Kruss tensiometer, according to the following steps. 

1 . A buffer solution is prepared at the predetermined pH. A biomacromolecule stock solution 
is prepared by dissolving the biomacromolecule in the buffer solution. An additive stock 
solution is prepared by dissolving the additive in the buffer solution. 

2. In FIG. 5, various solutions 120 are prepared at the predetermined temperature, usually 
room temperature, having the additive concentration 157 and different biomacromolecule 
concentrations 156 starting from 0 and increasing as long as the solutions 120 maintain a clear 
appearance. The solutions 120 are kept at the predetermined temperature for a few hours. The 
surface tension or surface pressure 164 of each solution 120 is measured by putting it into the 
Kruss tensiometer until the surface tension or pressure reaches a constant value. The chamber 
of the tensiometer is saturated with pure water vapor to maintain a uniform humidity, and the 
solutions 1 20 are kept at the predetermined temperature. 

The surface tension or surface pressure 164 is recorded and plotted against the 
biomacromolecule concentration 156. The critical magnitude 165 of the biomacromolecule 
concentration 156 occurring at the critical point 162 is registered. It corresponds to the surface 
tension or surface pressure 164, undergoing the transition 162 from the changing response 161 
to the substantially constant response 163 as the biomacromolecule concentration 156 
increases. The critical point 162 determined in this way is a solubility value of the 
biomacromolecule at the additive concentration 157. 
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Step 2 is repeated to obtain solubility points at different additive concentrations 157. 
Steps 1 and 2 are repeated to find the solubility values at different additive concentration 1 57 
and different pH. 

Second Embodiment of the Invention: Aggregation Boundary Condition 

In the second embodiment of the invention, the measurement of the aggregation 
boundary condition is carried out, as illustrated in FIGs. 4 and 5, by the following steps. 
L The buffer solution is prepared at the predetermined pH. The biomacromolecule stock 
solution is prepared by dissolving the biomacromolecule in the buffer solution. The additive 
stock solution is prepared by dissolving the additive in the buffer solution. 
2. Various solutions 120 are prepared having the biomacromolecule concentration 156 and 
different additive concentrations 157. As already stated, the biomacromolecule solution 120 in 
FIG. 4, used in the second embodiment for obtaining the aggregation boundary condition, 
should be fairly dilute. In the particular experiments described here the protein concentrations 
130 are e.g. 0.01, 0.1 or 1.0 mg/ml. Immediately after mixing each of the solutions 120, the 
surface tension is measured with respect to time and recorded. The measurement is stopped 
when the surface tension in each solution reaches a constant value. The obtained data are 
analyzed according to the steps outlined in FIGs. 7—15 described below. 
Step 2 is repeated for different pH. 

In the particular application, and without limitation to the scope of the invention, the 
experiment was carried out at the room temperature of 23 C. In FIGs. 4 and 5 the 
biomacromolecule 130 concentration is taken 1 mg/ml of lysozyme, the pH quantity is 4.5 and 
the buffering agent is 50 mM sodium acetate. In FIGs. 7 — 10 the additive is the salt NaCl, and 
the variable quantity 158, taken to be the additive concentration 157, is varied from 0 to 2.4 M 
in steps of 0.4 M. In FIGs. 7 and 8 the assembly parameter 164 is taken to be the surface 
tension s. In FIGs. 9 and 10 the assembly parameter 164 is the surface pressure p = s — so, 
where so is the surface tension of the solvent in the absence of biomacromolecules. For clarity, 
and without limitation to the scope of the invention, in FIGs. 7 — 10 only three values of the 
additive concentration 157 are used to illustrate the second embodiment of the invention: 0.0 
M NaCl 1 82, 0.4 M NaCl 1 84 and 1 .6 M NaCI 1 86. 

FIG. 7 shows plots 191, 192, 193 of the assembly parameter 164 taken to be the surface 
tension s, against time 180, for, respectively, the three additive concentrations 182, 184 and 
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186, taking the additive to be/the salt sodium chloride (NaCl) . We see that the corresponding 
assembly parameters 164, taken to be the surface tensions, remain unchanged when 
equilibrium times 194, 195 and 196 are reached. Each equilibrium time f eq corresponds to 
an equilibrium surface tension s cq . 

FIG. 8 shows plots 197, 198, 199 of the assembly parameter 164 against the square 
root of the time 180 for, respectively, the three additive concentrations 182, 184, 186. FIG. 9 
shows plots 202, 204, 206 of the logarithm of the time derivative of the assembly parameter 
(In dp/dt) 164, taken to be the surface pressure p> against the surface pressure p, for, 
respectively, the three additive concentrations 182, 184, 186. FIG. 10 illustrates time- 
dependence profile 220 of the macromolecule behavior against time 1 80, which expression in 
the present application is given by ln(l- p/p Gq ), where In is the natural logarithm, p is the 
surface pressure and p^q is an equilibrium surface pressure. The equilibrium surface pressure is 
given by p^ = s eq - so , where Seq is the equilibrium surface tension. 

In plots 212, 214, 216 of the configurational rearrangement expression 220 the 
equilibrium surface pressures p^ are equal to the surface pressure evaluated, respectively, at 
the equilibrium times t eq 194, 195, 196 in FIG. 7. 

We notice that each plot in FIGs. 8, 9 and 10 exhibits three line segments identifying 
the steps 142 144 and 146 in FIGs. 4 and 5: a first straight line segment 221 identifies the 
diffusion step 142, a second straight line segment 222 identifies the penetration step 144, and a 
third straight line segment 223 identifies the rearrangement step 146. In FIG. 10 the diffusion 
time fdiff is calculated by the inverse of the slope of the first straight line segment 221, the 
penetration time is calculated by the inverse of the slope of the second straight line 
segment 222, and the rearrangement time t m is calculated by the inverse of the slope of the 
third straight line segment 223, so that the integration time is obtained as a sum of the last two, 

Ant ~ fpen ^arr- 

FIG. 1 1 shows a plot 224 of the diffusion time f diff , a plot 226 of the penetration time 
'pen, a plot 228 of the rearrangement time and a plot 230 of the integration time t int = f pcn + 
fair, against the additive concentration 157. In this experiment the additive concentration 157 is 
taken in a range from 0 to 2.4 M NaCl in steps of 0.4 M NaCl. A histogram in FIG. 12 
illustrates how the crystallization coefficients 232 relate to the aggregation boundary 
condition. The crystallization coefficients 232 = /int/tam are calculated for all pairs of 
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obtained values for the integration time tint and the diffusion time f diff , and are shown against 
the additive concentration 157 in the above range. 

Crystallization Window 

Subsequently crystallization experiments were carried out in order to quantify the 
aggregation boundary condition in terms of the crystallization coefficient k^t 232 in FIG. 12. 
The experimental conditions used included the aforementioned range of the additive 
concentration 157 (FIGs. 1 1 and 12) as well as the remaining applied experimental conditions 
referring, e.g. to pH and temperature. The biomacromolecule concentration used for the 
purpose of crystallization was allowed to vary. The crystallization results are classified in three 
domains 234 listed below and included in FIG. 12, where the crystallization coefficient 232 is 
shown against the additive concentration 157. 

Domain A: The crystallization coefficient 232 is below approximately 4.0. Neither 
crystallization nor aggregation occurs in the entire range of the biomacromolecule 
concentrations used for crystallization, at the applied experimental conditions, e.g. additive 
concentration 157. 

Domain B: The crystallization coefficient feiyst 232 is above approximately 4.0 and below 
approximately 8.5. Crystallization occurs for certain biomacromolecule concentrations, at the 
applied experimental conditions, e.g. the additive concentration 157. 

Domain C: The crystallization coefficient Abryst 232 is above approximately 8.5. Only 
aggregation, but no crystallization, occurs at the applied experimental conditions, e.g. additive 
concentration 157, regardless of the biomacromolecule concentrations used for crystallization. 

Thus the aggregation indicator 232 employed to define the aggregation boundary 
condition in FIG. 12 lies below approximately 8.5, and the crystallization window is 
approximately given by the range 4.5 < < 8.5. An empirical result of our experiments is 
that in practice, a diffusion time that is short enough to enable crystallizing particles to 
approach one another and group together is reflected in a minimum value of = 4 or 4.5. 
That means, that below about feryst = 4 or 4.5 crystallization is usually not observed. The value 
of the crystallization coefficient indicating the onset of amorphous aggregation is the upper 
bound, empirically found to lie around Abiyst = 8.5 or 9. 

« 

Example 1. Experimental Result from Solution of the Protein Lvsozvme in Sodium Chloride 
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A combination of the crystal equilibrium condition in FIG. 1A or FIG. 6 and the 
aggregation boundary condition in FIG. IB or FIG. 12 are used to illustrate how crystallization 
and aggregation conditions can be predicted. As regards the crystal equilibrium condition, the 
point 171 on the solubility curve 170 of FIG. 1A or FIG. 6 corresponds to a pair of the additive 
(salt) concentration 157 that equals 1 M NaCl), and the biomacromolecule (lysozyme) 
concentration 1 56 (that equals 4 mg/ml). The crystal equilibrium condition is prescribing that 
crystallization can occur when the salt and lysozyme concentrations are taken above this pair 
of values. Thus the equilibrium lysozyme concentration equals 4 mg/ml when the additive 
concentration 157 is 1 M NaCl. The crystallization lysozyme concentration must therefore 
exceed 4 mg/ml. 

As regards the aggregation boundary condition, FIG. IB or FIG. 12 shows that when 
the salt concentration 157 is taken above around 1.9 M NaCl, no crystallization may be 
expected to occur, regardless of the lysozyme concentration used for crystallization; in other 
words any lysozyme concentration used for crystallization in this case has a very high 
probability of resulting in the amorphous aggregation. The crystallization coefficient feryst 232 
associated with the pair: salt concentration 157 equal to 1 M NaCl, lysozyme concentration 
156 equal to 1 mg/ml in FIG. IB or FIG. 12 is approximately 4.8. Therefore, when the salt 
concentration is 1 M NaCl, it should be possible for the crystallization to take place. That 
means that it should be possible to find some range of lysozyme concentration to use for the 
purpose of successful crystallization. Indeed, the crystallization does occur when the lysozyme 
concentration used for crystallization in this case is above 20 mg/ml. 

Example 2. Experimental Result from Solution of Protein Lysozyme in Ammonium Sulfate 

FIG. 13 shows a plot 240 of the assembly parameter 164 taken to be the surface tension 
s 9 against time 1 80, when the biomacromolecule solution 242 includes the protein lysozyme 
and the additive salt ammonium sulfate (NEUhSOi*. The biomacromolecule concentration 156 
is 1 mg/ml, and the additive concentration 157 is 2 M. This experiment resulted in the 
following values for the diffusion time, the penetration time, the rearrangement time, and the 
integration time: t d m = 833.33 sec, f pen = 2320 sec, fair = 5537 sec, f int = 7857 sec. The 
corresponding value of the crystallization coefficient 232 in FIG. 12 was 9.43. According 
to the aggregation boundary condition established in the present method, crystallization 
experiments using the aforementioned lysozyme and ammonium sulfate concentrations should 
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lead to the amorphous aggregation 128 in FIG. 3. Reports in the scientific literature confirm 
our result that lysozyme cannot be crystallized from the ammonium sulfate solution. However 
use of the method employing the second virial coefficient B22 leads to the incorrect prediction 
that the crystallization 126 in FIG. 2 of lysozyme from ammonium sulfate should be possible. 

Example 3. Experimental Result from Solution of Protein Concanavalin A in Ammonium 
Sulfate 

FIG. 14 shows a plot 244 of the assembly parameter 164 taken to be the surface tension 
s 9 against time 180, when the biomacromolecule solution 246 includes the protein 
concanavalin A (about 23 kDalton) and the additive salt ammonium sulfate (NEU^SO^ The 
biomacromolecule concentration 156 is 0.1 mg/ml and the additive concentration 157 is 1 M. 
This experiment resulted in the following values for the diffusion time, the penetration time, 
the rearrangement time, and the integration time: tdis = 277.78 sec, f pen = 1436.78 sec, = 
675.67 sec, /j n t = 2112.45 sec. The corresponding value of the crystallization coefficient 232 
kcryst in FIG. 12 was 7.605. According to the aggregation boundary condition established in the 
present method, crystallization experiments using the aforementioned concanavalin A and 
ammonium sulfate concentrations should lead to the crystallization 126 in FIG. 2. Reports in 
the scientific literature confirm our result that concanavalin A does indeed crystallize from the 
ammonium sulfate solution 246. No predictions based on the method employing the second 
virial coefficient ,822 are available in the literature. 

Example 4. Experimental Result from Solution of Protein Bovine Serum Albumin in 
Ammonium Sulfate 

FIG. 15 shows plots 251, 252, 253, of the assembly parameter 164 taken to be the 
surface tension s, against the time 180. The biomacromolecule solution 250 includes the 
protein Bovine Serum Albumin (BSA, 66 kDalton) and the additive salt ammonium sulfate 
(NHU)2S04. The biomacromolecule concentrations 156 are, respectively, 0.01 mg/ml, 0.1 
mg/ml 254 and 1.0 mg/ml, and the additive concentration 157 is 52% saturated ammonium 
sulfate. This experiment resulted in the following corresponding values of the crystallization 
coefficient 232 in FIG. 12: 24.6, 24.8 and 23.8. According to the aggregation boundary 
condition established in the present method, crystallization experiments using the 
aforementioned BSA and ammonium sulfate concentrations should lead to the amorphous 
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aggregation 128 in FIG. 3. Reports in the scientific literature however indicate that 
crystallization does indeed occur. This discrepancy is attributed to the tendency of BSA to 
dissociate and unfold on the surface 152 of the solution 120 in FIG. 4, i.e. in this case the 
solution/air interface. As a consequence of this tendency, the surface tension measurements 
prescribed in the present method offer an unreliable picture of the BSA crystallization 
window. 

Conclusion 

The prediction of the crystallization window through the embodiments of the present 
invention features the crystal equilibrium condition, with reference to FIGs. 1 A and 6, and the 
aggregation boundary condition, with reference to FIGs. IB and FIG. 12. Both embodiments 
are applicable for biomacromolecules that are hot prone to severe unfolding at the surface 152 
of the solution 120 (FIG. 4). The crystal equilibrium condition follows from the 
biomacromolecule solubility curve 170 (FIG. 6), and defines the lower bound of the 
biomacromolecule concentration to be used for crystallization. 

The boundary aggregation condition refers to the solvent experimental parameters. It is 
applicable when the biomacromolecule solutions used to determine the crystallization 
coefficient 232 (FIG. 12) are dilute. The results are summarized for various experiments 
in Table 1, where the employed biomacromolecules are proteins. The first column lists the 
protein concentrations 156 in FIGs. 7 — 15 used for predicting the crystallization window by 
determining the aggregation boundary condition. The protein concentrations used for the 
crystallization in the case of Iysozyme are allowed to vary in a range 5—100 mg/ml (not 
reported in Table 1). The second column lists the salt concentration 157 (FIGs. 7-15) used for 
determining the aggregation boundary condition, and also featuring in the crystallization 
experiment. The third column lists the pH, and the fourth column lists the resulting 
crystallization coefficient 232 (FIG. 12). The fifth column lists the outcome of the 
crystallization experiments in terms of the domains A, B or C, as predicted by the aggregation 
indicator of the present method. These predictions are compared with the predictions obtained 
from the method used in the prior art employing the second virial coefficient .822- Both 
methods are compared with experimental results in the last two columns of Table 1, whenever 
such results are available. 

Table 1 
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Protein, Salt, pH Crystallization B22 Surface tension 

concentration concentration result, prediction prediction 

(mg/ml) (M) domain (Fig. 12) 



lysozyme, 
1.0 


sodium chloride, 
<0.3 


4.5 


0.8 


solution, A 


correct 


correct 


lysozyme, 
1 0 


sodium chloride, 
0 3 2 0 


4.5 


4—7 


crystal, B 


correct 


correct 


lysozyme, 
1.0 


sodium chloride, 
>2.0 


4.5 


8—12 


aggregation, C 


correct 


correct 


lysozyme, 
1.0 


ammonium sulfate, 
2.0 


4.5 


9.5 


aggregation, C 


incorrect 


correct 


concalavine A, 
0.1 


ammonium sulfate, 
1.0 


7.0 


7.6 


crystal, B 


unavailable 


correct 


BSA, 

0.01,0.1, 1.0 


ammonium sulfate, 
52% saturated 


62 


24 


aggregation, C 


unavailable 


incorrect 



We see from Table 1 that the aggregation boundary condition in FIG. 12 prescribes that 
when the crystallization coefficient 232 is larger than approximately 8.5, the 
biomacromolecule has a high probability to undergo the amorphous aggregation 128 in FIG. 3 
at the employed experimental conditions involving, e.g., specifications of solvent, additive, 
their concentrations or the pH. Hence the aggregation indicator defines successftdly a 
relatively narrow crystallization window 4.5 < Abryst < 8.5, or a broader crystallization window 
4.0 < Acyst < 9.0, for biomacromolecules, with the following provisions: first that the 
biomacromolecules are able to maintain their conformations at the solution/air interface; and, 
second, relatively low biomacromolecule concentrations, e.g. not exceeding approximately 1 
mg/ml, are used in the prescribed experiments to determine the aggregation boundary 
condition. 
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1. A method for predicting a crystal equilibrium condition for biomacromolecule 

crystallization and for crystallizing a biomacromolecule, comprising 

setting up at least one biomacromolecule solubility experiment comprising the steps of 

a) preparing a solution of the biomacromolecule in a solvent, the solution having a 
biomacromolecule concentration, 

b) selecting a variable quantity, 

c) selecting an assembly parameter, 

d) monitoring a response of the assembly parameter while varying the variable quantity in a 
suitable way so that the response exhibits a transition, 

e) obtaining an equilibrium biomacromolecule concentration based on the transition, 

f) defining a crystal equilibrium condition according to which a biomacromolecule 
crystallization concentration exceeds the equilibrium biomacromolecule concentration, 

and crystallizing the biomacromolecule. 

2. The method as claimed in Claim 1, wherein the solution has further a pH and a temperature, 
and the variable quantity is one of the biomacromolecule concentration, the pH and the 
temperature. 

3. The method as claimed in Claim 2, wherein the solution further comprises an additive, the 
solution has an additive concentration, and the variable quantity is one of the 
biomacromolecule concentration, the pH, the temperature and the additive concentration. 

4. The method as claimed in Claim 1 , wherein the solution has a surface. 

5. The method as claimed in Claim 4, wherein the biomacromolecule is not prone to unfolding 
at the surface of the solution. 

6. The method as claimed in Claim 1, wherein the assembly parameter lis one of a density, a 
conductivity, a detergency and an osmotic pressure. 

7. The method as claimed in Claim 4, wherein the assembly parameter is one of a surface 
tension and a surface pressure. 

8. The method as claimed in Claim 2 or Claim 3, wherein the transition is associated with a 
critical magnitude of the variable quantity. 

9. The method as claimed in Claim 2 or Claim 3, wherein the transition is between a changing 
response of the assembly parameter and a substantially unchanging response of the assembly 
parameter. 
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10. The method as claimed in Claim 2 or Claim 3, wherein the transition is associated with a 
critical magnitude of the variable quantity, and further wherein the transition is between a 
changing response of the assembly parameter and a substantially unchanging response of the 
assembly parameter. 

11. The method as claimed in Claim 10, wherein the substantially unchanging response 
corresponds to a substantially minimal value of the assembly parameter. 

12. The method as claimed in Claim 10, further defining the crystal equilibrium condition in 
terms of the critical magnitude, wherein the crystal equilibrium condition prescribes that no 
crystallization occurs when the variable quantity is smaller than the critical magnitude. 

13. The method as claimed in Claim 12 wherein the variable quantity is the biomacromolecule 
concentration, and consequently the equilibrium biomacromolecule concentration equals the 
critical magnitude. 

14. The method as claimed in Claim 12 wherein the variable quantity is not the 
biomacromolecule concentration, and consequently the equilibrium biomacromolecule 
concentration equals the biomacromolecule concentration. 

15. The method as claimed in Claim 1, wherein the biomacromolecule to be crystallize^ is a 
protein. 

16. The method as claimed in Claim 15, wherein the protein has a weight less than 200 
kDalton. 

17. The method as claimed in Claim 16, wherein the protein is one of a lysozyme and a 
concanavalin A. 

18. The method as claimed in Claim 1, wherein the biomacromolecule to be crystallized is a 
polypeptide. 

1 9. The method as claimed in Claim 1 , wherein the biomacromolecule to be crystallized is a 
nucleic acid. 

20. The method as claimed in Claim 1, wherein the biomacromolecule to be crystallized is a 
virus. 

21. The method as claimed in Claim 1, wherein the biomacromolecule to be crystallized is a 
virus fragment. 

22. The method as claimed in Claim 3, wherein the additive is a salt. 

23. The method as claimed in Claim 3, wherein the additive comprises organic molecules. 

24. The method as claimed in Claim 3, wherein the additive comprises polymers. 
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25. A method for predicting a crystal equilibrium condition for protein crystallization and for 
crystallizing a protein, comprising 

setting up at least one biomacromolecule solubility experiment, comprising the steps of 

a) preparing a solution of the protein in a solvent, the solution further comprising an 
additive, the solution having a protein concentration, an additive concentration, a pH and a 
temperature, the solution having a surface, the surface having a surface tension and a surface 
pressure, the protein being not prone to unfolding at the surface, 

b) defining an assembly parameter to be one of the surface tension and the surface pressure, 

c) selecting a first variable quantity and a second variable quantity from the protein 
concentration, the additive concentration, the pH and the temperature, 

d) varying the first variable quantity in a suitable way so that the assembly parameter 
exhibits a transition between a changing response and a substantially unchanging response, 
wherein the substantially unchanging response corresponds to a first substantially minimal 
value of the assembly parameter, the transition being associated with a first critical 
magnitude of the first variable quantity, 

e) varying the second variable quantity in a suitable way so that the assembly parameter 
exhibits a transition between a changing response and a substantially unchanging response, 
wherein the substantially unchanging response corresponds to a second substantially minimal 
value of the assembly parameter, the transition being associated with a second critical 
magnitude of the second variable quantity, 

f) constructing a solubility curve comprising points, each point being a pair of the first 
critical magnitude and the second critical magnitude, in order to assist in defining a crystal 
equilibrium condition, 

g) obtaining an equilibrium protein concentration and defining the crystal equilibrium 
condition which is based on the solubility curve, and which prescribes that crystallization 
occurs when the first variable quantity exceeds the first critical magnitude of the pair, and the 
second variable quantity exceeds the second critical magnitude of the pair, 

and crystallizing the protein using a protein crystallization concentration exceeding the 
equilibrium protein concentration. 

26. The method as claimed in Claim 25, where in step (c) the protein concentration is one of 
the first variable quantity and the second variable quantity, and hence in step (g) the 
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equilibrium protein concentration is correspondingly one of the first critical magnitude and the 
second critical magnitude. 

27. The method as claimed in Claim 25, where in step (c) the protein concentration is not one 
of the first variable quantity and the second variable quantity, and hence in step (g) the 
equilibrium protein concentration is the protein concentration. 

28. The method as claimed in Claim 25, wherein the protein is one of the lysozyme and the 
concanavalin A and the additive is a salt 

29. A method for predicting an aggregation boundary condition for biomacromolecule 
crystallization and for crystallizing a biomacromolecule, comprising 

setting up at least one aggregation boundary condition experiment comprising 

a) preparing a solution of the biomacromolecule, 

b) selecting a variable quantity, 

c) selecting an assembly parameter, 

d) measuring the assembly parameter at different times, 

e) registering an equilibrium a5sembly parameter 

f) deriving a crystallization coefficient from the equilibrium assembly parameter, the 
crystallization coefficient being associated with the variable quantity, 

g) using an aggregation indicator to define an aggregation boundary condition for the 
biomacromolecule, the aggregation boundary condition prescribing that an aggregation 
occurs when the crystallization coefficient associated with the variable quantity is larger than 
the aggregation indicator, 

and crystallizing the biomacromolecule. 

30. A method for predicting an aggregation boundary condition for biomacromolecule 

crystallization and for crystallizing a biomacromolecule, comprising 

setting up at least one aggregation boundary condition experiment comprising 

a) preparing a solution of the biomacromolecule in a solvent, the solution having a 
biomacromolecule concentration and a surface, the surface having a surface pressure, 

b) selecting a variable quantity, 

c) obtaining the surface pressure at different times, while varying the variable quantity, 

d) recording a time dependent equilibrium surface pressure which is associated with the 
variable quantity, 
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e) formulating a time-dependence profile based on the equilibrium surface pressure, which 
is associated with the variable quantity, 

f) deriving from the time-dependence profile a crystallization coefficient of the 
biomacromolecule, that is associated with the variable quantity, 

g) obtaining from the crystallization coefficient an aggregation indicator in order to define 
an aggregation boundary condition for the biomacromolecule, the aggregation boundary 
condition prescribing that an aggregation occurs when the crystallization coefficient 
associated with the variable quantity is larger than the aggregation indicator, 

and crystallizing the biomacromolecule. 

31. The method as claimed in Claim 30, wherein the biomacromolecule is not prone to 
unfolding at the surface of the solution. 

32. The method as claimed in Claim 30, wherein the solution further has pH and a 
temperature. 

33. The method as claimed in Claim 30, wherein the biomacromolecule concentration is in the 
range 0.01 - 1.2 mg/ml. 

34. The method as claimed in Claim 30, wherein the solution further comprises an additive 
and the solution has an additive concentration. 

■ 

35. The method as claimed in Claim 32, wherein the variable quantity is one of the 
biomacromolecule concentration, the pH and the temperature. 

36. The method as claimed in Claim 34, wherein the variable quantity is one of the 
biomacromolecule concentration, the additive concentration, the pH and the temperature. 

37. The method as claimed in Claim 30, wherein the step of deriving the crystallization 
coefficient comprises the steps of 

obtaining a diffusion time of the biomacromolecule, 
obtaining an integration time of the biomacromolecule, 

dividing the integation time by the diffusion time to obtain the crystallization coefficient of 
the biomacromolecule, that is associated with the variable quantity. 

38. The method as claimed in Claim 30 wherein the time-dependence profile is given by ln(l- 
p/peq), where In is the natural logarithm, p is the surface pressure and is an equilibrium 
surface pressure. 

39. The method as claimed in Claim 38, where the step of deriving the crystallization 
coefficient comprises the steps of 
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constructing a plot of the time-dependence profile against a time, 

identifying on the plot of the time-dependence profile a first substantially straight linear 
segment, a second substantially straight linear segment and a third substantially straight linear 
segment, where the second substantially straight linear segment is later in the time than the 
first substantially straight linear segment and the second substantially straight linear segment is 
later in the time than the third substantially straight linear segment, 
equating a diffusion time to an inverse slope of the first substantially straight linear segment, 
equating a penetration time to an inverse slope of the second substantially straight linear 
segment, 

equating a rearrangement time to an inverse slope of the third substantially straight linear 
segment, 

adding the penetration time and the rearrangement time to obtain an integration time 
dividing the integration time by the diffusion time to obtain the crystallization coefficient of 
the biomacromolecule, that is associated with the variable quantity. 

40. The method as claimed in Claim 30, wherein the biomacromolecule to be crystallized is a 
protein. 

41. The method as claimed in Claim 40, wherein the protein has a weight less than 200 
kDalton. 

42. The method as claimed in Claim 41, wherein the protein is one of a lysozyme and a 
concanavalin A. 

43. The method as claimed in Claim 30, wherein the biomacromolecule to be crystallized is a 
polypeptide. 

44. The method as claimed in Claim 30, wherein the biomacromolecule to be crystallized is a 
nucleic acid. 

45. The method as claimed in Claim 30, wherein the biomacromolecule to be crystallized is a 
virus. 

■ 

46. The method as claimed in Claim 30, wherein the biomacromolecule to be crystallized is a 
virus fragment. 

47. The method as claimed in Claim 34, wherein the additive is a salt. 

48. The method as claimed in Claim 34, wherein the additive comprises organic molecules. 

49. The method as claimed in Claim 34, wherein the additive comprises polymers. 

50. The method as claimed in Claim 30, wherein the aggregation indicator is below 9. 
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51. The method as claimed in Claim 30, wherein the aggregation indicator is below 8.5. 

52. The method as claimed in Claim 30, wherein the aggregation indicator is in a range from 4 
to 9. 

53. The method as claimed in Claim 30, wherein the aggregation indicator is in a range from 
4.5 to 8.5. 

54. A method for predicting an aggregation boundary condition for protein crystallization and 
for crystallizing a protein, comprising 

setting up at least one aggregation boundary condition experiment comprising 

a) preparing a solution of the protein in a solvent, a salt, and a suitable buffer, the 
solution having a salt concentration, a protein concentration in a range 0.01 — 1.2 mg/ml, a 
pH and a temperature, the solution having a surface, the surface having a surface pressure, 
the protein not being prone to unfolding at the surface of the solution, 

b) obtaining the surface pressure at different times, while varying the salt concentration, 

c) recording a time-dependent equilibrium surface pressure, which corresponds with an 
equilibrium time, and which is associated with the salt concentration, 

d) formulating a time-dependence profile, which is given by ln(l-p//?cq), where In is the 
natural logarithm, p is the surface pressure and p eq is an equilibrium surface pressure, and 
which is associated with the salt concentration, 

e) constructing a plot of the time-dependence profile against a time, 

f) identifying on the plot a first substantially straight linear segment, a second 
substantially straight linear segment and a third substantially straight linear segment, where 
the second substantially straight linear segment is later in the time than the first substantially 
straight linear segment, and the third substantially straight linear segment is later in time than 
the second substantially straight linear segment, 

g) equating a diffusion time to an inverse slope of the first substantially straight linear 
segment, 

h) equating a penetration time to an inverse slope of the second substantially straight 
linear segment, 

i) equating a rearrangement time to an inverse slope of the third substantially straight 
linear segment, 

j) adding the penetration time and the rearrangement time to obtain an integration time 
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k) dividing the integration time by the diffusion time to obtain the crystallization 
coefficient of the protein, that is associated with the salt concentration, 

g) obtaining from the crystallization coefficient an aggregation indicator in order to define 
an aggregation boundary condition for the protein, the aggregation boundary condition 
prescribing that an aggregation occurs when the crystallization coefficient associated with the 
salt concentration is larger than the aggregation indicator, the aggregation indicator being in a 
range from 4.5 to 8.5. 

55. The method as claimed in Claim 54, wherein the protein is one of a lysozyme and a 
concanavalin A. 

56. The biomacromolecule crystallized according to any one of the Claims 1 — 24 and 30 — 53. 

57. The protein crystallized according to any one of the Claims 25 — 28, 54 and 55. 
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